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Abstract 

Multi-lepton signatures appear in many new physics searches at the Large Hadron 
Collider. We here consider WWZ production with subsequent leptonic decay of the 
three vector bosons as a SM source of multi-lepton events. We have calculated the 
next-to-leading order QCD corrections for the full pp — > 6 lepton production cross 
sections in hadronic collisions. Results have been implemented in the form of a flexible 
parton-level Monte-Carlo program which allows to calculate the QCD corrections for 
arbitrary distributions and acceptance cuts. 



1 Introduction 

With the start of the CERN Large Hadron Collider (LHC) a trove of new data is antici- 
pated, which can be used to search for new physics and also to further probe the Standard 
Model (SM). For the interpretation of these data, precise predictions for both the desired 
signal processes and backgrounds are needed, and this necessitates the calculation of next- 
to-leading order (NLO) QCD corrections. In order to determine cross sections for non-trivial 
acceptance cuts, it is most useful to cast these NLO calculations into fully flexible Monte- 
Carlo programs which can calculate cross sections as well as distributions. 

WWZ production with subsequent leptonic decay of the vector bosons is a background 
to supersymmetric processes with several leptons in the final state. In addition, it is an 
excellent probe of the quartic electroweak W + W~jj, W + W~Z~f and W + W~ ZZ couplings. 
Constraints for some of these couplings are already available from the Large Electron Positron 
collider (LEP) at CERN [I] and analogous measurements have been suggested for proton 
antiproton collisions at the Fermilab Tevatron collider [2J. However, the LHC will be able 
to improve these measurements considerably and therefore more accurate predictions are 
needed [3]. As we will see, QCD corrections increase the WWZ cross section by more than 
70% and, thus, any quantitative measurement of quartic couplings will have to take QCD 
corrections into account. This sizable enhancement is not surprising in view of the similarity 
with diboson production [U [5]. Furthermore QCD corrections to ZZZ production at the 
LHC have already been calculated and increase the LO cross section by about 50% [6]. 

We here present first results on the NLO QCD corrections to the full 2^6 process 
pp — > u e e + fi~ t~ t + (or any other combination of leptons from three distinct families). 
All resonant and non-resonant matrix elements as well as the spin correlations of the final 
state leptons are included in our calculation. For simplicity, we neglect any identical lepton 
effects which might appear when using the results for final states with leptons from one or 
two generations only. The calculation is performed with the Catani-Seymour subtraction 
algorithm [7] and uses the virtual amplitudes derived in Refs. [H [9]. Further details of the 
calculation and the checks we performed will be given in section 2. First determinations of 
cross sections, distributions and .fT-factors are presented in section 3. 

2 The Calculation 

We have considered in our calculation the full set of Feynman graphs for the process 
pp — > i/fj £\ ve 2 £3 £3 up to order a s a 6 . This includes the Higgs contribution and all off- 
shell diagrams. However, interference terms due to identical particles in the final state have 
been neglected. Including such effects at LO we find that this is an excellent approximation: 
LO cross sections change by less than 0.1% when interference terms are integrated over the 
Breit Wigner peaks. Interference terms show strong cancellations between contributions 
below and above the Breit Wigner peaks, but even their absolute values contribute at the 
few percent level only, which is below the scale variation of our final NLO cross sections. 
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Figure 1: Some representative tree-level Feynman diagrams of the process pp — > 4£ + 2 v. 
They show the three different topologies appearing in this calculation. 

We also neglect any fermion mass effects. In particular, the Higgs Yukawa coupling to tau 
leptons is set to zero. Because of the large number of Feynman diagrams we use the helicity 
method of Ref. [10] for the evaluation of the matrix elements. 

Some representative tree- level Feynman diagrams are given in Fig. [TJ The 181 tree- 
level graphs can be grouped into three distinct topologies. In Fig. dh three vector bosons 
with subsequent decay are emitted from the quark line. The polarization vectors of these 
vector bosons are the decay currents describing the respective decay leptons, and these decay 
currents appear in many different Feynman graphs. In order to speed up the calculation, 
they are determined numerically at the beginning of the evaluation of the matrix elements 
for a given phase space point and reused wherever they appear. In Fig. [Tb two vector bosons 
are attached to a quark line and then decay into two or four leptons. All the Feynman 
graphs for a 4-lepton decay can again be combined to an effective polarization vector. For 
all subprocesses like uu — > 4£ + 2u or dd — > A£ + 2u these polarization vectors are the 
same for one specific phase space point and they do not depend on the quark polarization. 
Furthermore they appear in various Feynman diagrams. In our code these polarization 
vectors are therefore calculated once per phase space point, stored and reused wherever 
possible. The last topology is the one where only one vector boson is attached to the quark 
line. The polarization vector corresponding to the "decay" of this virtual vector boson can 
be calculated once per phase space point, stored and reused. The method of precalculating 
effective polarization vectors renders our code for the Born processes about 4 times faster 
than a direct evaluation with MadGraph jl lj . 

The full NLO cross section consists of real emission contributions, virtual contributions 
and a finite remainder from the factorization of collinear singularities into the parton dis- 
tribution functions (pdfs). The real emission contributions have to be integrated over an 
(m+l)-particle phase space with m = 6 leptons, and the other two contributions only involve 
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an m-particle phase space. 
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(2.1) 



The virtual and the real emission contribution are separately infrared divergent and only 
their sum gives a well defined finite result. We have used the dipole subtraction algorithm 
proposed by Catani and Seymour [7] in order to handle these divergencies in our Monte-Carlo 
program. 

The real emission matrix elements can be obtained from the Born level matrix elements 
by either attaching a gluon to a quark line or having a gluon in the initial state which then 
splits into a quark anti-quark pair. The method of precalculating the effective polarization 
vectors for leptonic decays, described above for the tree-level diagrams, has also been used 
for the more complicated real emission diagrams. Here an increase in computational speed 
by about a factor of 12 is reached. 

The "virtual" contribution consists of the square of tree-level diagrams and the inter- 
ference between tree-level diagrams and the virtual one-loop diagrams. In the calculation 
of the one-loop diagrams, three different types of contributions corresponding to the three 
topologies in Fig. [T] appear. For the simplest topology, with one vector boson attached to 
the quark line, only vertex corrections appear. In this case the one-loop contribution is 
proportional to the Born matrix element. The second topology, Fig. [lb, leads to propagator 
corrections, vertex corrections and boxes. In the calculation, the sum of these corrections to 
one specific tree-level subamplitude is grouped together and will be called a boxline contri- 
bution in the following. The most challenging topology, Fig. [TJi, leads to quark propagator 
corrections, vertex corrections, boxes and pentagons. As in the previous case the sum of all 
these corrections to one specific tree level Feynman graph is grouped together and called 
pentline contribution in the following. 

The boxline contribution has the same structure as in the vector boson fusion process 
qq — > Vqq, which was considered in Ref. [8]. Similarly the pentline contribution is obtained 
by crossing from the results of Ref. [H], where QCD corrections to the vector boson fusion 
processes qq — > qqVV were determined. The singular contributions of all of these three 
types of virtual contributions are proportional to the Born matrix element and the complete 
virtual one-loop contribution for the three topologies is given by 
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where Mb denotes the full Born amplitude, s is the square of the partonic center of mass 
energy, i.e. it corresponds to the invariant mass of the 6-lepton system, and My is the 
finite part of the virtual boxline and pentline amplitudes, which are obtained by crossing 
and analytic continuation from the results of Refs. [El E] . In the calculation of the boxline 
routine the usual Passarino-Veltman tensor decomposition [12j is applied, whereas in the 
pentline routine the method proposed by Denner and Dittmaier [13] has been used for the 
tensor coefficients of pentagon diagrams. 
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Since the pentagon routines are quite time consuming we have employed a method sug- 
gested in Ref. [9] to reduce the magnitude of the true pentagon contribution. It is possible 
to shift the polarization vectors or respectively the decay currents of the vector bosons by 
terms proportional to their momenta, 



The pentagons contracted with the momenta instead of the polarization vector (terms pro- 
portional to xy) can then be expressed in terms of boxes via Ward identities for the loop 
integrals and the magnitude of the remaining pentagon contribution is reduced to a numer- 
ically less challenging level. In practice we choose 



which means that the shifted polarization vectors have zero time component in the center 
of mass system of the W-pair. 

We have checked our calculation at numerous levels. All the matrix elements for the 
LO process and the real emission part have been compared with MadGraph [IT] output. For 
individual matrix elements we find agreement at the 10 -15 level. In addition, the total LO 
cross section has been checked against MadEvent [IT] and HELAC [H]. For the real emission 
part the LO process pp — > 4£ + 2v + j has been compared against MadEvent. Results agree 
within the statistical accuracy of the Monte-Carlo runs (0.5% for MadEvent and 0.2% for 
HELAC). 

The finite collinear terms have been checked by exploiting the fact that they are generic 
for all Drell-Yan type processes and hence should be exactly the same for WW produc- 
tion. We have independently programmed NLO QCD corrections for WW production and 
compared our results with MCFM |5]. We find full agreement for different factorization and 
renormalization scales. As a further check we have taken advantage of the fact that we 
can integrate terms proportional to the Born matrix element either together with the real 
emission part, by integration over the (m+l)-particle phase space, or together with the 
virtual part, by integrating over the m-particle phase space. Resulting cross sections are 
independent of this choice of procedure. 

For the virtual contributions Ward identity tests have been implemented as a flag for 
numerically unstable evaluation of tensor coefficients. In the case of pentagon diagrams, for 
example, the pentline contribution can be expressed in terms of boxes if one of the external 
polarization vectors is replaced by its momentum. We discard the pentline contribution at 
phase space points where the two ways of calculating these terms disagree by more than 
10%. The resultant loss in cross section is negligible: extrapolating its size we estimate the 
resultant relative error on the NLO cross section to be below 10~ 4 . 

Many of the critical elements of our calculation also appear in the related process pp — > 
ZZZ for which NLO QCD corrections have recently been presented in Ref. [6], albeit without 
further decay of the leptons and neglecting any Higgs contributions. As a last and very 
powerful test we have modified our code to describe ZZZ instead of WWZ production and 
checked against the results of Ref. [5]. Using the same pdfs and a s (mz) = 0.119 at NLO we 
find complete agreement for all of the given renormalization and factorization scales. 



e'y = x v q v + e'y 



(2.3) 
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3 Results 



The complete next-to-leading order calculation has been implemented in the framework 
of a fully flexible Monte-Carlo program, VBFNLO, which then has been used to determine the 
QCD corrections to WWZ production at the LHC For the generation of our results we set 
the electroweak parameters to the following values: 

m w = 80.419 GeV m z = 91.188 GeV 

Gp = 1.16639 ■ 1(T 5 GeV" 2 m H = 120 GeV. (3.1) 

The other electroweak parameters, a -1 = 132.507 and sin 2 (9w) — 0.22225, are calculated in 
the program by using LO electroweak relations. The default value for the renormalization 
and factorization scale is the invariant WWZ mass, which is given by 

Hf = Vr = V (Ph + Pi 2 + Ph + Pu + Pvi + Pu 2 ) 2 - (3.2) 

For scale variation studies of the total cross section, we also consider fixed renormalization 
and factorization scales, which we set as //p = hr = £ • mz- We use the CTEQ6M parton 
distribution with a s (mz) = 0.118 at NLO and CTEQ6L1 for the LO calculation [15J. All 
fermions are treated as massless and we do not consider contributions involving bottom 
and top quarks. The CKM matrix is approximated by a unit matrix throughout, which is 
appropriate when neglecting fermion mass effects in a neutral Drell-Yan type process. 

In order to keep most of the cross section we apply only minimal cuts to the final state 
leptons. This means that besides a cut on the minimal transverse momentum and the 
maximal rapidity of the charged leptons we only require that the invariant dilepton mass, 
ran, of any combination of the charged leptons is larger than 15 GeV, thus steering clear of 
the virtual photon singularity in 7* — > £ + C,~ at low m«. Specifically, we require 

pr t > 10 GeV \yt\ < 2.5 m u > 15 GeV (3.3) 

in all subsequent figures. All results given here have been calculated for the process pp — > 
u e e + pT t~ t + , i.e. interference terms due to identical particles have been neglected. In 
order to obtain cross sections for the phenomenologically more interesting case of final states 
with four electrons and/or muons, we have multiplied these results by a combinatorial factor 
of 8 in all figures. 

In the left panel of Fig. [2] the scale dependence of the cross section at LO and NLO 
is shown. The LO scale dependence severely underestimates the uncertainties of the LO 
cross section due to the fact that no as appears in the LO calculation and that the pdfs 
are determined in a Feynman-x range of small factorization scale dependence. The right 
panel of Fig. [2] shows the scale dependence of the different contributions to the NLO cross 
section. The virtual contribution is split up into one part proportional to the Born matrix 
element, one part with boxline contributions and a third part which includes only the true 
pentline contribution after the reduction described in equations (12.31) . f)2.4p . This pentagon 
contribution amounts to 1 - 2% of the total NLO cross section. The real emission contribution 
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Figure 2: Left: Scale dependence of the total cross section for pp — > WWZ — > 4£ + ^ T at LO 
and NLO for m# = 120 GeV and the cuts of Eq. A3.3\) . The factorization and renormalization 
scales are taken at a fixed value which is varied in the range from 0.5 • m>z to 10 • rriz- Right: 
Same as in the left panel but for the different NLO contributions. 



in the plot includes the additional finite terms from the factorization of collinear singularities 
into the pdfs. The scale dependence of the NLO cross section is mainly due this real emission 
contribution. 

The .ff -factor, defined as the ratio of the NLO cross section over the LO cross section 
strongly depends on the factorization and renormalization scale. It is shown in Fig. [2] for a 
fixed scale choice, [Ir = fip = £ ■ rriz- One finds variations of the rT-factor between 1.7 (for 
large scales) and 2.2 (for small scales). A qualitatively similar scale dependence is found if 
instead of the Z-mass the invariant 6-lepton mass as given in Eq. (13.21) is taken as reference 
scale choice. The mean value of this dynamical reference scale is substantially larger than 
the Z-mass and the K-factor for /ip = /ir = mwwz is 1.74. Since the WWZ invariant mass 
is the scale of the Drell-Yan type subdiagrams in Fig. 1, it may be considered the most 
natural scale choice in the present calculation. We leave a more thorough investigation of 
scale dependence of the NLO cross section to a future publication [IB] , however. 

The size of the NLO QCD corrections shows a marked phase space dependence. Two 
examples are shown in the remaining figures. In the left panel of Fig. [3] the transverse 
momentum distribution of the highest-pr charged lepton is shown for the cuts given in 
Eq. (13.31) . The right panel shows the differential K-factor, defined as 

da^o/dx 

K = d*">/dx ■ ( 3 ' 4 ) 

The rT-factor increases with px by almost a factor 2, indicating that a simple multiplication 
with a constant overall i^-factor would seriously change the shape of lepton px distributions. 
A somewhat more benign behavior is found for the 4-lepton invariant mass distribution which 
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Figure 3: Left: Transverse momentum distribution of the highest-px charged lepton for 
rriH = 120 GeV, u>f = Hr = rnwwz as given in Eq. ( 13. 2|) and the cuts given in Eq. h3.3\) at 
LO and NLO. Right: Differential K-factor as dehned in Eq. A3.4\) for the two distributions 
in the left panel. 
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Figure 4: Left: Invariant mass distribution of the four charged leptons for ttlh = 120 GeV, 
= = rnwwz and the cuts given in Eq. ( 13.3)) at LO and NLO. Right: Differential 
K-factor as dehned in Eq. ( 13.4)) for the two distributions in the left panel. 
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is shown in Fig. HI Here the differential i^-factor varies between 1.4 and 1.9 and most of 
the rise is in the threshold region. For 4-lepton invariant masses above 250 GeV a constant 
i^-factor would be an adequate approximation for the distribution. 

4 Conclusions 

WWZ production with subsequent leptonic decays is an important source of multi-lepton 
events and of great interest for the measurement of quartic electroweak couplings at the LHC. 
In this letter we have presented first results on the NLO QCD corrections to the process 
pp — > vn x £\ £2 i>i 2 £3 For the WWZ invariant mass as renormalization and factorization 
scale, the overall .ff -factor is about 1.7 and therefore has to be taken into account in any 
analysis involving WWZ production. 

The scale dependence of the NLO cross section is substantially larger than the variation 
observed for the LO results. This can be quantified by increasing and lowering the renormal- 
ization and factorization scale by a factor of 2 around 3-mz as central value. At LO the scale 
dependence is very small with a variation of less than ± 1.5%, whereas at NLO variations 
of ± 5% appear. The NLO uncertainties which are indicated by the scale dependence are, 
thus, typical for a NLO QCD prediction, while the LO case must be considered as anoma- 
lous. Indeed, the WWZ cross section provides another example where the scale variation 
of a LO cross section does not give a good estimate for the corrections due to higher order 
effects. 

The NLO QCD corrections do not only change the normalization of the total WWZ 
cross section, they also lead to substantial shape changes of distributions: the differential 
.fT-factors show a sizable variation over phase space. Reliable modeling of the lepton distri- 
butions arising from WWZ production at the LHC thus requires the inclusion of NLO QCD 
corrections. These corrections are now available in the form of a flexible parton level Monte 
Carlo program, which will be incorporated into the publicly available VBFNLO package [T7] 
in the near future. 
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